Reconstruction of the photon distribution in a micromaser 
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We suggest an iterative, maximum-likelihood-based, method to reconstruct the photon number distribution of 
the steady state cavity field of a micromaser starting from the statistics of the atoms leaving the cavity after the 
interaction. The scheme is based on measuring the atomic populations of probe atoms for different interaction 
times and works effectively using a small number of atoms and a limited sampling of the interaction times. 
The method has been tested by numerically simulated experiments showing that it may be reliably used in 
any micromaser regime leading to high-fidelity reconstructions for single-peaked distributions as well as for 
double-peaked ones and for trapping states. 
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Introduction — The one-atom maser or micromaser is per- 
haps the most relevant example of open quantum system in 
cavity quantum electrodynamics (CQED) HI]. Since its first 
experimental realization ]j3] this system has allowed to inves- 
tigate many fundamental aspects in quantum optics. The mi- 
cromaser dynamics results from the interplay of a coherent 
interaction between a beam of two-level atoms and a reso- 
nant cavity mode in the microwave domain, as described by 
the Jaynes-Cummings (JC) model and the dissipative pro- 
cess due to the contact of the cavity with the environment. 
At the steady state the radiation field inside the high-Q cav- 
ity may show highly non-classical features, as for example 
sub-Poissonian photon statistics ^ or quantum collapses and 
revivals In addition, states characterized by a truncated 
photon number distribution the, the so-called trapping states 
(TS) of the cavity field may be generated. These states 
show up only at very low temperature and may be affected 
by collective atomic interactions A- Under suitable pump- 
ing conditions, photon distribution at the steady state may 
also show two coexistingmaxima, that is the signature of first 
order phase-transitions ■ Operating the micromaser under 
pulsed regime and trapping conditions the generation of Fock 
states has been also reported ]i9|]. A micromaser was imple- 
mented also on a two photon transition ifioll and more recently 
a microlaser was operated in the optical regime 111 1 1. 



A crucial aspect of the micromaser is that, in order to pre- 
serve high Q values of the cavity, the cavity field is not ac- 
cessible to direct measurements. As a consequence, any in- 
formation on its properties must be inferred from the atoms 
leaving the cavity after the interaction. In fact, the relation 
between the atomic statistics and the properties of the cav- 
ity field has been theoretically and numerically investigated 
I jJO also including the back-action due to the atomic measure- 
ments ifisll . On the other hand, to the best of our knowledge, 
no method has been suggested to reconstruct the whole photon 
distribution by exploiting the complete information carried by 
the atoms leaving the cavity. In earlier experiments ^ U the 
atomic statistics was obtained by counting the number of ex- 
cited (ground) atoms in a time interval longer than the cavity 
lifetime, and then this frequency was compared with the the- 
oretical expression (see below) for the experimental set of pa- 



rameters. As a matter of fact, the photon distribution was not 
reconstructed from the measurements. In experiments lead- 
ing to TS ]0] the steady state photon distribution is composed 
only by few terms allowing a simple fit of experimental data, 
whereas in the experiments to generate Fock states ]i2] the cav- 
ity field state \n) is prepared by a pulse of n pump atoms and 
only one probe atom is measured to obtain the atomic inver- 
sion that ideally involves only one Rabi frequency. In this 
case, the advantage to measure only one probe atom is that of 
avoiding the cavity field state reduction due to repeated atomic 
measurements. 

In this letter, we suggest a method to reconstruct the whole 
steady state photon distribution of the cavity field starting 
from measurements of the statistics of probe atoms. The ba- 
sic idea is that atoms leaving the cavity after different inter- 
action times are carrying the complete information about the 
cavity field itself. Indeed, the method is based on measuring 
the atomic statistics for different interaction times and then 
estimating the photon distribution using maximum-likelihood 
reconstruction. As we will see, the method is very effective 
in any operating regime of the micromaser and allows reli- 
able reconstructions for single-peaked distributions as well as 
for multi-peaked ones and for trapping states. Remarkably, 
the method works effectively starting from the statistics of a 
small number of atoms and a limited sampling of the inter- 
action times. As a consequence, the atoms used to probe the 
cavity field are only slightly perturbing the steady state, which 
itself depends on the interaction time of the pump atoms, i.e. 
the method can be used on-line with experiments. We also 
notice that at the steady state, the cavity field density matrix 
is diagonal in the Fock number basis, and thus the reconstruc- 
tion of the photon distribution corresponds to the full quantum 
state reconstruction. On the other hand, the characteristics 
of the micromaser spectrum 11411 are related to the decay of 
off-diagonal elements of the cavity field density matrix in the 
transient regime. 

Photon distribution at the steady state — A schematic dia- 
gram of the micromaser setup is given in Fig.[T]where a beam 
of two level atoms, excited in the upper Rydberg level of the 
maser transition, continuously and resonantly pump a high-Q 
microwave cavity mode. The cavity temperature is kept as low 
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as 0.5K in order to have a small number of thermal photons. 
The velocity of the atoms can be selected so that the interac- 
tion time tint between each atom and the cavity mode can be 
selected with high precision. The atomic flux has a Poisso- 
nian distribution with a mean pump rate R. The state of the 
atoms leaving the cavity can be detected by field ionization 
techniques. 
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FIG. 1 : Schematic diagram of the micromaser setup. 

The atomic decay rate and the cavity decay rate 7 are 
taken such that tint <C ^ ^ 7~^. Under the above 
conditions only one atom interacts with the cavity mode each 
time, thus realizing a perfect JC interaction, and a steady state 
regime for the cavity field can be obtained. If gp denotes the 
(diagonal) steady state density operator of the cavity field, the 
photon distribution p„ = (Oint, ^cx7 ?^tii) = {n\gF\n) can 
be expressed as [jst] 

" (N^^/m) sin^ (eint^/mJN^) + nth 

Pn=PoY[ Vt ' 

1 + nth 

where po is a normalization constant, Nex = R/j the ef- 
fective pump rate, and 8int = gtint V ^cx the dimensionless 
pump parameter, g being the atom-cavity coupling constant. 

A striking consequence of Eq. ^ is the existence of trap- 
ping states of the cavity field f^. In the limit of nth the 
distribution p„ vanishes at photon numbers nq (q = 1, 2, ...) 
such that 9int = qT^\/Ncx/ (1 + nq). The TS correspond to 
narrow dips which appear in the stationary mean photon num- 
ber (N) = ^Pn as a function of the pump parameter 8int- 
Another interesting form of p„ can be obtained if the pump 
parameter is set to Qiat = § corresponding to maximum am- 
plification (MA) regime of the micromaser In this case, p„ 
has a shape like that of a coherent state with the same mean 
photon number Finally, close to Qhit — Stt and multiples 
thereof, the photon distribution p„ assumes a double-peaked 

SP) structure corresponding to a first-order phase transition 
i. 

When the system is at steady state, the probability Pe to 
find one atom in the excited state after its interaction with the 
cavity field for a time tk is given by 
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n=0 
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1 + COS (nV" + 1) 



where Pfe = Pe(Tfe) and Tk = gtk, is the dimensionless inter- 
action time. Eq. (|2]i provides a link between the experimen- 
tally measurable statistics of the probe atoms and the (inac- 
cessible) photon distribution of the cavity field. 



Reconstruction of the photon distribution — At a first sight, 
Eq. Q seems to provide a scarce piece of information about 
the photon distribution p„ of the micromaser However, if 
the atomic statistics is recorded for a suitable set of values of 
the interaction times, then the information is enough to recon- 
struct the full photon distribution. As we will see, the inver- 
sion of Eq. Q, i.e. the reconstruction of p„, may be obtained 
by maximum-likelihood estimation upon a suitable truncation 
of the Hilbert space. 

The reconstruction scheme proceeds as follows: the micro- 
maser is pumped until it has reached the steady state for a 
fixed set of parameters. Then we stop the atomic pump flux 
and a probe atom prepared in the excited state is sent through 
the cavity in a time much shorter than the cavity photon life- 
time 7^^. The velocity of this probe atom may be adjusted 
in order to vary the interaction time in a given range. After 
the interaction with the (steady state) cavity field the probe 
atom is detected. We denote by fk = Mk/N^ the experi- 
mental frequency of probe atoms found in the excited state 
after an interaction time Tk, Nx being the total number of 
atoms sent through the cavity with interaction time Tk. Of 
course, since atom detection modifies the cavity field state, 
every probe atom is followed by pump atoms to restore the 
steady state field. In the following, we assume that the values 
of interaction times for the probe atoms Tk, k = 0, . . . , n,- 
are uniformly distributed between a minimum value tq and a 
maximum one , which, in turn, are determined by the max- 
imum and minimum velocities allowed by the specific exper- 
imental implementation. 

Eq. (|2]i is a statistical model for the parameters pn that can 
be solved by maximum-likelihood (ML) estimation. We as- 
sume that the photon distribution can be truncated at the n-th 
term (i.e. p„ is negligible for n > n) and, without loss of gen- 
erality, that Nx is independent on k. The global probability of 
the sample i.e the log-likelihood (with normalized P^) of the 
detected data reads as follows: 
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ML estimates of Pn are the values maximizing the log- 
likelihood L. Since the model is linear and the unknowns p„ 
are posi tive the solutio n can be obtained using an iterative pro- 

can be 



cedure Us , 3 12 , 18 ] . Indeed, the equations 
written as 
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1 Vn = 0, 



(3) 



Then, by multiplying both the sides of Eq. (O by p„, we get 
a map Tp„ = p„, whose fixed point can be obtained by the 
following iterative solution 



P: 



where pn''' is the value of p„ evaluated at the h-th iteration. 
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and Pi'')(rfc) =E „ Ckn Pn . Eq. (|4| is usually referred to as 
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the expectation-maximization solution of ML, and is known 
to converge unbiasedly to the ML solution. As a matter of 
fact, Eq. dUi provides a solution once the initial distribution 
Pn^ is chosen. In our simulated experiments we start from the 
uniform distribution = (1 + in the interval [0, n], 
though any other distribution such that X^n?'""' ~ 1' 
Pn^ 7^ Vrt, would be appropriate as well. Indeed, the initial 
distribution is slightly affecting only the convergence rate and 
not the precision at convergence 11911 . 
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FIG. 2: Reconstruction of the photon number distribution from 
Monte Carlo simulated experiments for different steady state micro- 
maser regimes. On the left we report the probability Pe{r) of find- 
ing an atom in the excited state as a function of the interaction time 
r (as obtained from the reconstructed distribution, solid line) com- 
pared with the actual frequencies observed in the simulated experi- 
ments (crosses). On the right we show the reconstructed photon dis- 
tribution (crosses) compared with the theoretical one (histograms). 
The micromaser parameters are A'^ox ~ 25.0, nth ~ 10~^ and (a) 
TS regime, Gint/Tr = 2.5; (b) MA regime, Oint/Tr = 0.5; (c) DP 
regime, Qmt/n — 2.18. In all the simulated experiments tq = 0.5, 
Ur = 200, and A^^ = 200, Nu = 1000. 

Monte Carlo simulated experiments — Reliability and ac- 
curacy of the present method have been tested by an exten- 
sive set of numerically simulated experiments, correspond- 



ing to different micromaser steady state regimes. As a fig- 
ure of merit to assess the accuracy of the reconstructed dis- 
tribution p„ , i.e. the similarity to the actual distribution 

Pn of Eq. ([Hi, we consider the fidelity G ~ ^„ \J pn^ Pn- 
In Fig. |2] we show the simulated experimental data for the 
measurement of Pe(T), generated by Monte Carlo technique, 
and the comparison between the theoretical photon distribu- 
tions and those obtained by ML estimation. We consider 
as interesting examples the TS, MA and DP steady state 
micromaser regimes. In these regimes the photon number 
distribution is sub-Poissonian, nearly Poissonian, and super- 
Poissonian, respectively. In order to better appreciate the ac- 
curacy of our reconstruction method we also report (see Ta- 
ble ^ the first two moments of the cavity field distribution, 
i.e. the mean photon number n = (a^a) and the Fano factor 
F = [((a^a)^) — (a^a)^] /{a^a), a being the mode operator 
of the cavity field and (• • • ) = Tr [g-p ■ ■ ■] denoting ensemble 
average. As it is apparent from TableUa very good agreement 
is obtained for all the considered regimes between the values 
obtained from the reconstructed distributions and the actual 
ones. 



TABLE I: Mean photon number and Fano factor of the reconstructed 
distributions of Fig.[2]compared with the theoretical values. 
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TS 


2.5 TT 


2.52 


0.22 


2.53 


0.21 


99.73 


MA 


0.5 TT 


24.38 


1.02 


24.34 


1.04 


99.94 


DP 


2.187r 


7.85 


4.05 


7.76 


4.03 


99.72 



Being our reconstruction method based on an iterative so- 
lution an important aspect to keep under control is its conver- 
gence. In Fig. [3] we show the fidelity of the reconstruction as 
a function of the number Ur of sampling interaction times and 
the number N^. of measures for each interaction time, for dif- 
ferent numbers of iterations A'^it. As it is apparent from Fig. [3] 
the fidelity increases with both Ur and Nx and it reaches an 
asymptotic value which actually depends on the choice of the 
other parameters. Of course, also the number of iterations A'^it 
affects the fidelity value at convergence. Notice, however, that 
the reconstruction is already very accurate with a number of 
iterations A^it = 100. It is worth noticing that the number 
of sampling times tIt- cannot be increased at will, since it is 
limited by experimental constraints. In order to check the sta- 
tistical reliability of the algorithm we report the results from 
repeated (simulated) experiments. The error bars in Fig. [3] 
are obtained by averaging over one hundred simulated exper- 
iments. 

A question arises on whether the present method could be 
effectively employed with a small number of atoms and a lim- 
ited sampling of the interaction times. This is of course a cru- 
cial aspect concerning its possible implementation in a realis- 
tic scenario. We found, by means of an extensive set of simu- 
lated experiments, that the answer is positive and that accurate 
reconstructions may be obtained using realistic values of the 
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interaction times. This features, together with its accuracy and 
fast convergence, make it suitable for being used on-line with 
experiments. The method has been tested by numerically sim- 
ulated experiments showing that it may be reliably used in any 
steady state regime of the micromaser leading to high-fidelity 
reconstructions for single-peaked distributions as well as for 
double-peaked ones and for trapping states. 
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FIG. 3: (a): Fidelity of the reconstruction G as a function of the 
number of sampling times at fixed number of data Nx = 200 for 
each time value, (b) Fidelity of the reconstruction G as a function of 
the number Nx of data for each time value at fixed number Ut — 200 
of sampling interaction times. Both plots refer to the case of TS state, 
i.e. the reconstruction reported in Fig. [2ji. In both plots we report 
the fidelity for different values of the number of iterations Nit, from 
bottom to top: Na = 100, 200 and 1000. The error bars are obtained 
by averaging over one hundred simulated experiments. 



parameters. In Figs. |4] we report, as an example, the results 
of a simulated experiment, corresponding to that of Fig.|2c), 
now performed with rir ~ 40, N^. — 30 and A'it — 50. As 
it is apparent from the plots, the reconstruction is still very 
accurate despite the fact that the total number of observations 
has been dramatically decreased. 
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FIG. 4: Reconstruction of the photon number distribution as in 
Fig.|2jc) but with reduced sampling parameters Nx = 70, n-r = 50, 
A'^it = 300, and n = 25. For these values of parameters we 
have n = 11.55, F = 2.34, n'"' = 11.41, F'"' = 2.35, and 
G = 98.14%. 

Summary and conclusions — We have suggested a novel 
iterative method to reconstruct the full photon distribution of 
the cavity field of a micromaser at the steady-state starting 
from the statistics of the probe atoms leaving the cavity af- 
ter different interaction times. Our methods works effectively 
using a small number of atoms and a limited sampling of the 
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